Static and Dynamic properties of a Fermi-gas of cooled atoms near a wide Feshbach 

resonance. 
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We introduce the Global Spin Model to study the static and dynamic properties of the ultracold 
fermionic gas near the broad Feshbach resonance. We show that the problem of molecular produc- 
tion, in a single-mode approximation, is reduced to the linear Landau-Zener problem for operators. 
The strong interaction leads to significant renormalization of the gap between adiabatic levels. In 
contrast to static problem the close vicinity of exact resonance does not play substantial role. Two 
main physical results of our theory is the high sensitivity of molecular production to the initial 
value of magnetic field and generation of a large BCS condensate distributed over a broad range of 
momenta in a wide range of parameters. We calculate the amplitude of the condensate as function 
of the initial detuning and the rate of the magnetic field sweep. 



I. INTRODUCTION 

In recent years there have been numerous achieve- 
ments in the area of ultra-cold atomic physics. The 
major experimental tool for it is the use of the Fesh- 
bach resonances (FR)ii^i^i^i^i^ii, which occurs when the 
energy of a quasibound molecular state becomes equal 
to the energy of two free alkali atoms. The magnetic- 
field dependence of the energy allows precise tuning 
of the atom-atom interaction strength in an ultracold 
gas^. Moreover, time-dependent magnetic fields can be 
used to reversibly convert atom pairs into weakly bound 
molecules^'^ i^^'^^i^^'^'^i^^'^^ . This technique has proved to 
be extremely effective in converting degenerate atomic 
gases of fermions^ i^'^i^^i^^i^^i^^i^^ and bosonsiiii^ii^ into 
bosonic dimer molecules. 

The main theme of this paper is the crossover from 
the BCS to BEC condensates. The weakly and strongly 
paired fermionic s-wave superfiuids were well understood 
by 60's, while the relation between the superfiuids were 
studied later by Eagles^^, Leggctt^, and Nozicres and 
Schmitt-Rink^^. These works treated the BCS state as 
a variational ground state and explicitly showed that 
there is no qualitative distinction between the BCS and 
BEC superfiuids. After the discovery of High- Tc mate- 
rials a possibility to describe the pseudogap phase as a 
long crossover between the BCS and BEC condensates 
renewed the interest to the phenomena— i^. It also be- 
come very important in the area of ultracold atomic gases 
where the transition/crossover can be directly observed 
in the molecular production experiments. 

Theoretical works on the molecular production can be 
roughly referred to two categories. The first is a phe- 
nomenology suggesting that pairs of molecules indepen- 
dently undergo Landau-Zener (LZ) transitionsi*^. The 
total number of molecules at the end of the process in 
these works is the LZ transition probability multiplied 
by the number of pairs. The most problematic issue 
in this approach is how to identify a pair in the gas 
of indistinguishable fermions. Direct calculation of the 
transition probability from a microscopic Hamiltonian, 



to the 4-th order in the interaction constani— , shows 
that, in contrast to the assumption of phenomenologi- 
cal works, the many-body effects are essential. Another 
category includes works based on a simplified model^^ 
in which molecules have one available state mimicking 
the condensat o^^i^^i^^ . Although numerical works in this 
category display a reasonable temperature dependence, 
they give no clear physical picture and detailed parame- 
ter dependence. 

In this article we consider a homogeneous gas in a box 
with the large volume V. The gas consists of Fermi iso- 
topes of alkali atoms such as ^Li, ^^Na, "'"K. The gas 
is subject to an external magnetic field such that the 
electronic Zceman energy strongly exceeds the tempera- 
ture, whereas the nuclear Zecman energy remains much 
smaller than the temperature. Therefore, the electronic 
spins of atoms are almost completely polarized along the 
magnetic field. The nuclear spins in the ground state are 
then aligned with respect to the electronic spin by the 
hyperfine interaction. In a typical experimen t ^^'^^ the 
magnetic field passes the Feshbach resonance. 

The origin of this resonance is the following. When two 
atoms with one electron in the outer shell each turn into 
a molecule, the two electrons form a spin singlet. This re- 
quires one of the electrons to fiip its spin, which costs the 
Zeeman energy for the electron. The resonance happens 
when Zeeman energy equals to the binding energy of the 
molecule. Neither exchange nor magnetic dipolar inter- 
action of electrons can produce such a change of the spin 
state. This process is mediated by the hyperfine interac- 
tion. The characteristic scale of the hyperfine interaction 
£/i / is in the range of niK and it is smaller or of the same 
order of magnitude as the molecular binding energy. It 
strongly exceeds the Fermi-energy of the gas which is in 
the range of 10^^ A'. Collisions in the s-wave channel 
of Fermi-atoms is allowed by the Pauli principle only if 
they have different states of nuclear spins. Therefore, ex- 
perimenters create an admixture of about equal numbers 
of the fermionic atoms with different nuclear states, for 
example atoms of ^'^K with the same total atomic spin 
9/2 but different spin projection quantum numbers —7/2 
and —9/2. Further we will describe this different states 
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by a pscudospin index a accepting two values "f , J, . 

The difference between the Zeeman energy and its res- 
onance value is called detuning e. It is the experimentally 
controlled parameter, which can be changed by changing 
external magnetic field. There are two different types 
of experimental procedures. In the first one, which wc 
call static, an equilibrium or mctastable state is estab- 
lished at a fixed value of e and the number of atoms 
and molecules as well as the distributions of their veloc- 
ities as function of detuning are studied. In the second 
one, which we call dynamic, the magnetic field sweeps 
passing the FR at which the detuning energy e is zero. 
After passing through the resonance a part of atoms are 
transformed into diatomic bosonic molecules. The ex- 
perimental and theoretical problem is the description of 
this process and finding the molecular production as a 
function of the magnetic field rate. 

A major difference between these experiments is in 
time scales. At the first, static, type of experiment the 
system is always at thermal equilibrium (quasiequilib- 
rium), while at the second, dynamic, type the gas is out 
of equilibrium - the equilibrating processes are too slow. 

A general feature of the equilibrium and dynamic 
state is the appearance of the condensates: the Bardeen- 
Cooper-Schrieffer (BCS) condensate on atomic side of the 
resonance and the Bose-Einstein (BE) condensate on its 
molecular side. The condensates have the same symme- 
try. Therefore, the static transformation from the BCS 
to BE condensate proceeds not as a sharp transition, but 
rather as an eventual crossover. 

The main subject of our paper is the BCS and BEC 
condensates densities and their correlations in both static 
and dynamic problems at condition of a broad resonance 
which will be shown to be equivalent to a strong interac- 
tion. As always the strong interaction problem is difficult 
to solve exactly. We propose a reasonable approximation 
which allows the analytical solutions of both, static and 
dynamic problems. 

Below we introduce the necessary notations and defini- 
tions. The dependence of the s-scattering length a on ex- 
ternal magnetic field B is determined by the well-known 
equation: 




a (B) = ao 1 



B-B^ 



(1) 



There are two different regimes for both static and 
dynamic transitions: narrow and broad Feshbach reso- 
nance. The Feshbach resonance is narrow when A is 
much smaller than the Fermi energy sp and broad in 
the opposite case. Thus, the resonance is broad if the 
dimcnsionlcss parameter 



(3) 



is large. In the opposite case the resonance is narrow. 
The static and dynamics of the narrow resonance was 
thoroughly analyzed theoretically by Gurarie and Radzi- 
hovsky to whose review^S we refer the reader. The case of 
the broad resonance was not analyzed since it is a strong 
interaction problem. Wc show that the condition of the 
broad resonance pI7|) allows simplification of the model 
despite of strong interaction between fermions. The key 
idea is a proper cut-off in the momentum space and ne- 
glect of the fcrmion dispersion. The resulting model is 
similar to the Dicke spin model for superradiance^. This 
model allows us to solve the static problem exactly. The 
complete spectrum and eigenstates are found. The solu- 
tion displays a crossover from BCS to BEC in a range 
of detuning near the FR. We find the density of the con- 
densates and their correlation as function of the detuning. 
We show that, contrary to the standard BCS solution, the 
scale of momenta in which the Cooper pair wave function 
is confined is (2mA) much larger than the Fermi en- 
ergy. Therefore, the size of the Cooper pair in real space 
is much less than the average distance between atoms. 
We discuss consequences of this fact for both statics and 
dynamics. 

The dynamic problem is reduced to the Landau-Zencr 
problem in the operator form which also can be solved 
exactly. One of conservation laws in our theory ensures 
that the increase of one of condensates proceeds only at 
expense of another. The main conclusions of the dynamic 
problem arc: i) Sensitivity of the molecular production to 
the initial value of the BCS condensate; ii) Conservation 
of a sum of weighted BCS and BEC condensate densities 
which indicates that one of condensates can grow only at 
the expense of another one. The dynamic problem was 
briefly described earlier in our article^. 



where oq is the scattering length far from the resonance; 
Bq is the resonance value of the field and B\ is the 
magnetic-field width of the resonance which appear in 
experimental measurements. An intrinsic energy scale 
generated by the field width is^S: 



fi2 



(2) 



We will show later that the value of A characterises the 
interaction in the Fermi-gas or the interaction of the BCS 
pairs with the BEC condensate. 



II. THE MODEL AND CONSERVATION LAWS 



We start with the Timmermansl Hamiltonian^: 
il = ^(e + £p)a^^ap^ + ^ cjq6]j6q -f 

-% (^«P + qT"-pi + ^qO-Pi^P+qt) (4) 

^ p,q 

Here aj,^^ are creation operators of fermionic atoms with 
momentum p, pscudospin cr, and kinetic energy £p = 
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p^/2m; 6]j are the creation operators of the bosonic di- 
atomic molecule with kinetic energy uiq = q^/4m; the 
position of the FR is controlled by the experimentally 
tunable detuning energy e, wich becomes time-dependent 
in the dynamic problem. The last term in ^ is the in- 
teraction with a coupling constant g. It describes the 
formation of a molecule from two atoms and the inverse 
process of dissociation into two atoms. 

As it was explained earlier, the coupling constant g is 
provided by the hyperfinc interaction. Its value can be 
estimated^ as g ~ ehf\/afj^, where en f is a character- 
istic hyperfine energy (about ImK) and a„i is the size 
of the diatomic molecule. The Hamiltonian ^ neglects 
nonresonant direct atom-atom and molecule-molecule in- 
teractions that near any FR arc subdominant to the res- 
onant scattering. To some extent the direct atom-atom 
interaction is taken into account by the value uq in cq. 
([!]) . We use the word "atoms" for both uncoupled atoms 
and atoms within molecules, and the word "fermions" for 
atoms that are not bound in molecules. Correspondingly, 
we denote the number of atoms as N and the number of 
fermions as Np = X]p a- cr'^p,'^- Since the diatomic 
molecule can be created by the interaction term in ^ 
only at expense of two fermions the Hamiltonian ([4]) con- 
serves the total number of atoms: 



N = Nf + 



(5) 



Therefore, in the Hamiltonian ([4]) the detuning energy 
can be transferred with the coefficient —2 to the energy 
of molecules. To find the connection between the cou- 
pling constant g in the Hamiltonian (jj]) and the value A 
defined by equation ([2]) we use the perturbation theory 
to the second order to eliminate bq and fej^ from the static 
Hamiltonian ([4]) and obtain the Hamiltonian for fermions 
only with 4-fermion interaction. Assuming that Sp and 
ujq can be neglected in comparison to e (this assumption 
will be justified later), the interaction Hamiltonian reads: 



Hi, 



9 

2Ve 



p.p',q 



'p+qT 



•J 



pX"-p'i'^p'+qT 



(6) 



Thus, the effective interaction constant for fermions is 
gp = g^/2e. It is negative (attraction) at negative e. The 
s-scattering length a is related to the interaction constant 
gp by a standard relation^ gp = A-Kh^a/m. Thus, the 
singular part of the scattering length is associated with 
the fermion-boson coupling constant g and the detuning 
energy e as follows: 



a(e) 



(7) 



Comparing equation ([7]) with the resonance term in 
equation ([T|) and identifying e ~ ^b{B — Bq), we ar- 
rive at a relation 2Biao^B = mg'^/ATrh?. Using ep = 
/i^(37r^n)^/'^/2m, where n is the density of fermions, and 



substituting this into equations (|2l3p we express the in- 
trinsic energy A and dimensionless parameter F in terms 
of the coupling constant g: 



m^g'^ 
167r2ft6 

9 9 

m'^g'^ 



1 



(8) 

^4„l/3 ^5/331/323/2 ' 

The Hamiltonian ([4]) is still too complicated. We show 
that the broad resonance condition 



F> 1 



(10) 



allows us to apply two simplifying assumptions. 

The first one of them is the single-mode approximation: 
we neglect all bosonic modes except the one with zero mo- 
mentum. The characteristic energy scale in statics and 
dynamics is A. The corresponding value of the scattering 
length following from equation ([5]) at such a detuning is 
a(A) ~ n-i/3/F, which gives the characteristic value of 
the gas parameter an^^^ ~ F~^ ^ 1. SmaU |e| ^ A/F 
correspond to the strong interaction regime in which the 
fermions are distributed in a broad range of momentum 

The maximal modulus of the scattering am- 
plitude / for such particles achieved in the unitary limit 
is / h/ V 2mA. The corresponding value of the gas 
parameter fn^/^ is again of the order F^^ ^ 1. At this 
conditions and at zero temperature almost all particles of 
the Bose-gas fall into the coherent condensate with zero 
momentum. We thus will substitute bq = ^o.q^o in the 
Hamiltonian ([¥]). This is the single-mode approximation. 

The applicability of the single-mode approximation to 
the dynamics seems less obvious since the formation of 
the condensate in the Bose gas requires some time which 
must be compared with the characteristic time of the 
fermion-boson transition. The kinetic of the condensate 
formation was theoretically studied by Kagan et al^t^ 
and by Gardiner et al^. In the article^ it was concluded 
that the formation of condensate starts after some cool- 
ing process which requires a delay time, typically few 
tens of ms in the laser cooled gases and then the con- 
densate grows rapidly. In the experiments^ even smaller 
delay time was found. The characteristic time of sweep- 
ing through the FR was of the same order of magnitude 
or longer. Besides, the experimenters started the sweep- 
ing in a close vicinity of the resonance in which at least 
the BCS condensate was sufficiently strong. As we will 
see later, in dynamics it produces the BEC condensate 
during the time characteristic for the Landau-Zener pro- 
cesses. 

The second approximation we apply is disregarding the 
fermionic dispersion in ^ in comparison to £. A typical 
value of bo is ^/N. In the broad resonance approximation 
g\/N/\/V = g\/n is much larger than ep ^ h^ii?/^ /m., 
gm^/n/h?!!?/^ ~ \/V > 1, see (|9]). With this motivation 
we neglect the kinetic energy of fermions in comparison 
to e. Though rather rough, this approximation retain es- 
sential features of the exact Hamiltonian. The character 
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and precision of this approximation will be discussed in 
more details later. 

These two approximations greatly simplify the Hamil- 
tonian ^ reducing it to the following form: 



Ps Pa 

bo OptaLpx +blY, apia-Pt 



(11) 



.9 

w 



The number of atoms conserved by the Hamiltonian pT 
reads: 



Ps 



2blbo 



(12) 



p,<T 



In Eqs. pip and we have introduced the cutoff mo- 
mentum Ps to account for the dropped fcrmionic disper- 
sion. The value of this momentum must be such that our 
approximation of neglecting the kinetic energy in com- 
parison to A is justified. Thus, a reasonable cut-off mo- 
mentum is: 



Ps = V 2mA 
Plugging eq. ([8|) into eq. p3|) . we find: 



Ps 



? ? 



(13) 



(14) 



The number of available states TV,,, which includes all 
the states with p < Ps, reads: 



Ns^V 



Ps 



(15) 



Eq. p4p gives the following result for the density of the 
available states — Ns/V: 



(16) 



29/2 . 37^5^12 

With precision of the coefhcient 1.09 the value 

A = 5V^ (17) 

coincides with A. The condition of the broad resonance 
is equivalent to the strong inequality: 



(18) 



Let us consider the whole Hilbert space as a set of pairs 
of conjugated states |p, t) and |— p, |). Each pair can be 
either empty, or single occupied, or doubly occupied. A 
single occupied pair is not changed by the Hamiltonian 
pT|) . The whole Hilbert space can then be split on two in- 
variant subspaces: the set of the single occupied pairs and 
the rest, where no state is single occupied. The Hamil- 
tonian is diagonal in the first subspace, so it induces no 



transitions in it. We thus consider only an invariant sub- 
space of the Hilbert space with all pairs either empty or 
doubly occupied. 

Besides of the occupation number of fermionic pairs, 
the state is fixed with the number of bosonic molecules 
iY„i. The Hilbert subspace invariant under the action of 
the Hamiltonian pT|) at a fixed total value of atoms N 
consists of states with given number Af,„ of molecules in 
the condensate and the rest of atoms occupying pairs of 
fermionic states. The number of occupied pairs is M = 
N jl ~ Nm- Summing the number of states with different 
possible M, we find for the dimensionality of the Hilbert 
space 'D{N,Ns) with given N and Ns'. 



N/2 

V{N,Ns)^ J2 



Nsl 



M=0 



M\iNs - My. 



(19) 



III. 



THE GLOBAL SPIN MODEL AND ITS 
HILBERT SPACE 



Following Anderson^S, we introduce spin operators: 

Spz = ^ (flpT^PT + "-pi"-pi - l) (20) 



(21) 



In the double-occupied or empty fermionic pairs subspace 
they obey the standard spin-1/2 commutation relations: 

[sp2,Sp'±] = ±(5pp'Sp±; [sp+,Sp/_] = 2^pp/Sp^ (22) 

so that the double occupied pair corresponds to Spz = +5 
and empty pair corresponds to Spz = ^5- Note that a 
single-occupied fermionic pair corresponds to a singlet 
spin state: all three spin operators (|201 BT|l turn such a 
pair to zero. 

The neglect of kinetic energy allows us to rewrite the 
Hamiltonian pip in terms of only global spin operators: 



Hs = 2eSz + ( boS+ + blS^ 



(23) 



(we have omitted a constant originated from the term — 1 
in equation (|20p ). where the global spin components 



'pz J 



5± = E' 



■p± 



(24) 



obey the standard commutation relationships: [S'z,S'±] = 
±S±; [S+S^] = 2Sz- We wiU caU the model with the 
Hamiltonian ^ the Global Spin Model (GSM). 

A subtle assumption incorporated in the derivation of 
the Hamiltonians pip . ([^5]) . and in the permutation re- 
lations for the global spin components is that the sum- 
mations in the two definitions (|24p run over the same 
range of momenta. It is not obvious since the summa- 
tion in the sum for Sz is limited by the condition Sp < |e|, 
whereas the summation in the second sum originated 
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from the boson-fcrmion interaction part of the Hamil- 
tonian is naturally cut off by the range of interaction. It 
means that g is a function of momentum vanishing at 
sufficiently large values of the momentum modulus. In 
the dynamical problem the Anderson spins rotate with 
the frequencies (e + ep) /h. Therefore they rotate coher- 
ently and enhance remarkably their effective field exerted 
to the moleculare amplitude only if Ep < |e|. We will 
see that such a coherence indeed takes place in the dy- 
namic problem. The contributions from larger values of 
momenta to S± are incoherent and mutually compensate 
their effect. Thus, for dynamic problem, the summation 
in the same momentum region in the two equations (|24[) 
is justified. 

It is not so obvious for the static problem. We will show 
that the static GSM has many qualitative features resem- 
bling what is expected for the initial Timmermannsi et 
al. model and observed in the experiment, though no- 
body solved the latter model exactly. In particular the 
GSM displays a crossover from BCS to BEG condensate 
with a large gap due to strong interaction in a broad 
vicinity of the Feshbach resonance. However, in quanti- 
tative details they are different. Besides, the GSM ([23]) 
has an additional symmetry and an additional conserved 
value which the Timmermans Hamiltonian does not pos- 
sess. Indeed the Hamiltonian ([25]) conserves not only the 
value of 



Sz + ^^0^0, 



equivalent to the number of atoms ((12 
spin S, where: 

S{S + l) = Sl + \ {S+S_ + S-S- 



(25) 



but also the total 



According to ([181) . '^z is always large and negative. The 
total spin S cannot be smaller than {S^l, it also cannot be 
larger than Ns/2. Thus, the allowed values of the total 
spin S are {N^ -N)/2<S< Ns/2. 

The number Af{S, Ng) of different multiplets with a 
fixed value of the total spin S which appear at the addi- 
tion of Ns spins 1/2 reads^ 



M{S,Ns) = 



Nsl{2S+l) 



(29) 



Each of these representations contains generally 25+1 
states, but only S — of them are allowed by the 

inequalities (j28p . The number of states Af{S,Sz) in the 



Hilbert space of the GSM with fixed values S and Sz 
reads: 



U{s,Sz) = e{Sz + s)e 



2 



N 
1 



Sz\M{S), (30) 



(26) 



where B(x) is the Heaviside step function. 

The Hamiltonian (P5)) formally coincides with the fa- 
mous Dicke Hamiltonian for the so-called superradiance 
problem^ij but it acts in a different Hilbert space. 



IV. SPECTRA AND EIGENSTATES OF THE 
GSM 

The Hamiltonian (|23p commutes with S and N . There- 
fore, its stationary states |^) at fixed values A^s, A^, and 
S can be represented by a superposition of states with 
fixed number of molecules Nm or S, =■ — A',,, : 



In the thermodynamic limit of large system it is possible 
to neglect 1 in comparison to S and consider 5*+ and S- 
as commutative values. The conservation laws (^51 1^5)) 
imply that the quantity of one of the condensates (BCS or 
BEG) can be increased only at the expense of the other. 

The cut-off in the momentum space define a finite- 
dimensional Hilbert space of states i?jv. What we call 
the GSM is this Hilbert space and the Hamiltonian 
acting in it. 

Each state in i?jv is a vector from a direct prod- 
uct of Ns spin 1/2 representations corresponding to 
the Anderson spins Sp. The combinatorial coefficient 
NsljMMNs - M)\ entering eq. (dH) can be treated as 
the number of possible distributions of M spins up and 
Ns — M spins down, i.e. the number of different states 
with the spin projection: 



(27) 



Since < M < N/2, the projection 5*2 at a fixed S is 
limited by inequalities: 



S<Sz<{N~ Ns)/2. 



(28) 



(31) 



whose amplitudes ^n„ 
equation: 



obey the stationary Schrodinger 



= -2€Nn'^N^ + 

-^T V^™ (S - Sz) {S + Sz + l)*iv„-i + (32) 

vV 



^ ViNm + 1){S + Sz) {S-Sz + l)*Ar„. + l 

This system is still rather complicated, but it is strongly 
simplified by the broad resonance condition (fT5)) . Indeed, 
due to this inequality S ~ —Sz « Ns/2, and it is possible 
to replace S — Sz and S — Sz + I in eq. ([55)1 by A^^. 
We should be more careful with the expression S + Sz 
since the two terms almost completely cancel each other. 
Employing the notation defined earlier A = g^/n^ (fT7|l 
and introducing new variables 



Sz 



-m: 



S N 
2+- 



m 



Nm-s, (33) 
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we arrive at a simplified version of equations (|33p : 



-2em*™ + Av/(s - m)(s + m + 1)4-^+1 + 
Av/(s + m)(s-m + l)1',„+i + 2se)«'™,(34) 

Since A'^m ^ and S* + > 0, the quantum number 
m runs from —s to s and < s < Tlie equation 

(p4)) is easily recognizable as generated by a reduced spin 
Hamiltonian: 



-2es, + 2As, 



(35) 



where Sx-, Sz are spin operators corresponding to the re- 
duced total spin s. This is a Hamiltonian of a spin s in 
the magnetic field 2V + tilted in the a-z plane at 
the angle 6* = — tan~^(A/e) to the z axis. 

The energy levels are labeled by two integers s and m 
(we have introduced tilde to distinguish projection to the 
axis tilted by the angle 9 to z— axis from projection to 
the z— axis): 



Es,n (e) = 2TOyPTA2 - 2se 



The spectrum (|36p possesses a symmetry: 

Esrh (e) ~ —Es-m (" e) 



(36) 



(37) 



Levels with the same s and different rh do not cross, but 
the levels with different s and m cross each other. The 
crossing of the levels (s, m) and (s', m') happens at a 
point: 



sign 



> TO 



\ m — ra J 



(38) 



1 



Besides of crossings each level (s, rh) at any e is A/'(5, Sz) 
fold degenerate as it is determined by eq. ([29|) . In terms 
of reduced spin variables it reads: 



N{N, s, rh) « e{s - m)e{s + in) 



27r(iV/2-2s)! 



(39) 

For a fixed s the state with minimal energy is (s, — s). 
The ground state corresponds to the maximal possible 
value s = N/4: and rh = —s. Its energy is: 



Er 



E 



N 



N/4, N/4: 



(40) 



The energy of the ground state is separated from the rest 
of spectrum by a finite energy gap S, which is determined 
by the following equation: 



(5 = 2 



V?+A2 + e6 



(41) 



The first excited state at negative e is the state with 
N 1 „„A ™ _ K ^ I- positive e it is 

N 



s = — — 1 and TO 
the state with s 



N/A and rh 



1. In the 



ground state the maximal spin s = N/A is oriented 
in the xz— plane at the angle 9 = tan^^(A/e) to the 
z— axis. This state can be thought of as a set of N/2 
spins 1/2, all directed along the same line. The cor- 
responding wave function is = riilfi^ ^IV"*); where 
= (cos0/2, sin0/2)J'. The operator Sz is represented 

by the direct sum: Sz = 'Yl!i=i Therefore, the aver- 
age value of TO in the ground state is equal to {m)Q — 

WszW = E^JiiHsim = fcose = Ne/AV^^T^. 
Employing the third equation we find the average 
number of molecules in the ground state: 



{N, 



N 



m/G 



+ A^ 



(42) 



It smoothly varies from at e = ~oo to N/2 at e = +oo. 
The width of the transition is determined by A. The 
value {Nm)Q — is an odd function of the detuning 
energy. The average number of fermions (Np) can be 
found from the conservation law Np + 2Nm = N . 

The fluctuation of the number of molecules in the 
ground state also can be calculated exactly. In- 
deed, ((AiV„)2)G = (V'|(s,)2|^) - {',},\szW using 

(V|(s.)2|^) =7V/8 + E5^'(V',|4lV'.)(V',|silV',> = A^/8 + 
cos^ {9)N{N - 2)/16. Thus, at large N, the quadratic 
fluctuation of the number of molecules reads: 



{W-^)^)g 



N A^ 
8 e2+, 



(43) 



The fluctuation ((AA^,„)2) is maximal at e = — at the 
FR, and is an even function of the detuning e. 

An important value is the BCS condensate amplitude 
{S+)g in the ground state. Using the first equation ([55]) 
and (IMl), we find that \{S+)g? ~ S^-{S^)g « {Ns/2)^~ 
{NJ2 - iV/4)2 - N,{sz)g » iNNs{l - cos 61) and thus 



1(5. 



+>g\ 



NN, 



sin 6*72 1 



(44) 



The ground state is symmetric with respect to any per- 
mutation of the Ns pair states, which is the Hamil- 
tonian symmetry operation. Therefore, all amplitudes 
('^pT'^-pi^ are equal and each of them is equal to 
yA/V/]Vs sin(0/2). A simple physical picture behind this 
distribution is the following. The total number of occu- 
pied pairs is equal to N sm^{9/2). They are equally dis- 
tributed among Ns available states. Therefore, the prob- 
ability to find a pair at any fixed state is N siT?{9/2)/Ns 
and amphtude to find this state is yjN/Ng sin(0/2) x e*'^p . 
This result coincides with indicated above if all ampli- 
tudes are coherent and have the same phase. According 
to this picture only the total number of pairs depend on 
detuning, but their distribution in momentum space re- 
mains unchanged: they are uniformly distributed over 
momenta from zero to Ps- It means that the size of the 
Cooper pair is always h/ps- it is much smaller than the 
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distance between gas particles n~'^/^ . The pairs arc com- 
pact, but still the collective interaction of the fermions 
inside the pairs play important role. The implicit as- 
sumption of the theory is that the size of the Bosonic 
molecule is much less than h/ps- Another important dif- 
ference between the Cooper pair and molecule is that in 
the former electron spins of atoms are parallel being po- 
larized by magnetic field, whereas in the latter the elec- 
tron spins are antiparallel. The broad distribution of the 
Cooper wave function in the momentum space is a conse- 
quence of the corresponding broad range of the hyperfine 
interaction. This property will persist in a more accurate 
theory which does not accept our simplifying approxima- 
tions. 

We notice from ([42|) and (gi]) that the value 



NsN,n + \{S+)\^ =NsN/2 



(45) 



is a detuning independent constant. 

Finally, we find the BEC-BCS correlation function 
(6o'S'+). We demonstrated earlier (see the derivation of 
eq. from (|33p ) that in the broad resonance approx- 
imation Ng ^ N the product of operators boS+ can be 
replaced by 



{boS+) = 



/K,s+. Thus, 



4 Ve^ + A2 



(46) 



This correlator vanishes at e = ±oo and has maximum at 
e = 0. An interesting relationship between the BEC-BCS 
correlator and the fluctuation of the number of molecules 
follows from cqns. (gS]) and (ge]): NsN{{ANm)'^)G = 
2{boS+f. 

The GSM displays BCS-BEC crossover in the range 
|e| ~ A near the Feshbach resonance. In this range 
the BCS condensate amplitude grows to the value ~ 
y/NsN/2 > N/2. This enhancement of the conden- 
sate is due to the distribution of the Cooper pairs over 
a wide range of momenta strongly exceeding the Fermi 
sphere. It indicates that the famous BCS exponentially 
small condensate docs not appear even at very large de- 
tuning exceeding A. In the BCS theory the condensate 
or energy gap is exponentially small not only due to the 
weakness of interaction, but also because the attraction 
range in the momentum space is very narrow. The latter 
condition is violated not only in the GSM, but also in the 
initial Timmermanns et al. model. 



V. DYNAMIC ATOMS-TO MOLECULES 
TRANSFORMATION 

In this section we consider the transformation of Fermi- 
atoms into molecules in the cold Fermi gas under the 
sweep of the magnetic field. In the Hamiltonian pTjl the 
parameter e is driven by magnetic field: e = ^b{H — Hq). 
At magnetic field sweeping, the detuning energy e de- 
pends on time. The time-dependent Hamiltonian regu- 
lating the process of atoms-to-molcculcs transformation 
reads: 



H = 2e{t)S, + 



bs+ + y^s^ 



(47) 



It formally coincides with the GSM Hamiltonian ([23]) 
with the only difference that the detuning energy e de- 
pends on time. Even the time-dependent Hamiltonian 
(|47)) commutes both with the operator Q given by (|25p . 
equivalent to the operator of the total number of atoms 
iV, and with the square of the total spin operator S"^ ([^5]) 
which both remain the integrals of motion. The Heisen- 
bcrg equations of motion are: 



-igS-\ hS- = -2ie{t)S- + 2igb'^S^ (48) 



hb 



where g = g/vV and dots denote the time derivatives. 
Generally these equations are non-linear. However, in 
the broad-resonance approximation Sz = —Ns/2, they 
become linear. Eliminating we arrive at an ordinary 
linear differential equation for the operator b: 



h^b + 2ihe{t)b - 



A^b 







(49) 



where A is defined in pT|) . 

Equation (g^ becomes the parabolic cylinder equation 
if e(t) is a linear function of time. In the LZ theory it 
describes the evolution of the probability amplitude to 
find the system in one of its two states. The value A in 
the corresponding LZ problem is the matrix element of 
transition between the two crossing terms, the so-called 
LZ gap. 

As we already argued, during the magnetic field sweep 
the Anderson spins with p <^ ps rotate coherently with 
almost the same angular velocity and enhance the effec- 
tive field acting on the BEC condensate. The spins with 
p > Ps rotate with substantially different angular veloci- 
ties and mutually cancel their contribution to the effec- 
tive field. Our model coarsens these features, neglecting 
the decoherence of spins with p < ps and completely ig- 
noring spins with p > p^. 

An important conclusion is that strong interaction 
renormalizes the LZ gap. The energy scale which ap- 
pears in perturbation theory is A^*^^ = gy^rS^. For a 
broad resonance A = g^/nl is much larger than A^'^^ and 
does not depend on the atomic density. 

Employing equations (|48|) . the general solution of the 
ordinary differential equation (|49p reads (further we put 
h=\): 

b[t) = u (i, to) 6 (to) - i~gv (t, to) S- (to) , (50) 
igS-it) = -ii (t, to) b (to) + igv {t, to) (to) , (51) 

where u(t,to) and v{t,to) are standard solutions of 
the same equation (|49|) satisfying the initial conditions 
u {to, to) = 1, ii {to, to) = and v {to, to) = 0, v {to, to) 
1. These solutions have the following properties: 



\uUA-'\u'\ 



^A^Id^H 



\uf+A''\vf 



u*v + A'^u*v — 0; 



uv 



1; (52) 

if*^e(t')dt' ^ (53) 
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The solution (|50l5ip completely determines the evolution 
of the number of molecules Nm{t) = {wb){t), the BCS 
condensate amplitude F{t) defined by equation F'^{t) = 
{S+S-){t), and the BCS-BEC coherence factor C{t) = 
{wS-){t). It is convenient to introduce the intensive val- 
ues: the molecular (or EEC) density nm{t) = Nm{t)/V, 
the BCS condensate density f{t) = g^jF{t)/V and the 
BEC-BCS correlator density c{t) = -~igC{t)/V^/'^. Then 
equations of motion read 

n„{t) = |M|2n,„o + |«|'/o + 2Re {u*vco) (54) 
f{t) = |itpn,„o + |i;p/o + 2Re(ri>co) (55) 
c{t) = uu*nmo + vv* /q + u*vco + uv*Cq, (56) 

where we also introduced n,„o = nmXto), /o = /(^o)j a-nd 
Co = c(io) ^ the initial values of the density of molecules, 
the BSC condensate density, and the BEC-BCS correla- 
tor density. Using ((52|) and summing eqs. ([54)1 and ([55)) . 
we find: 



+ = const. 



(57) 



which is a consequence of the conservation laws. Since 
for any state F'^{t) > 0, if there are no molecules in the 
initial state, their number Nm{t) can not exceed the value 
F^{to)/Ns at any time. Although equation ([57)1 has the 
same meaning as equation (|45|) it is more general as it 
shows that one of condensates can be extended only at 
the expense of another even at a finite sweep rate. 

Below we consider two experimentally most relevant 
situations: i) only fcrmions and no molecules; ii) only 
molecules and no fermions in the initial state. In both 
these cases the initial value C = {b^S-){to) = 0. In the 
case of no molecules in the initial state, so Nm{to) = 0, 
the general equations (|54l56p simplify to 

Nra{t) = -g^lvl^F^to); i~gC{t) = -~g''vv* F^{to) (58) 



The evolution of F{t) in this case is determined by ([57]) 
and ([55]) . Note that the coherence factor C{t) does not 
remain zero. 

In order to produce a reasonable fraction of molecules 
it is necessary to have a large condensate amplitude in 
the initial state with the only exception for the equilib- 
rium initial state and adiabatically slow sweep of mag- 
netic field. Due to the finite gap in the spectrum our 
model predicts that in this situation the system will adi- 
abatically follow the ground state. However, the initial 
state of the gas in the experiment is reached by the pump- 
ing of the ac electromagnetic field and it is plausible that 
it is not in equilibrium. 

A strong dependence of the final molecular production 
on the initial state (in particular on the value of the ini- 
tial magnetic field) explains why different experimenters 
obtain different fractions of molecules in the final state 
even in the adiabatic regim o^^i^^i^^i^^ . 

In experiments, which achieved a significant molecu- 
lar production, the initial state was indeed close to the 



FR, whereas the final state was rather far from the reso- 
nance. Thus, in a realistic experimental setup the initial 
value of e is small |eo| < ^ A and then e increases 
linearly with time. In this case one can put t^ = 0, 
and e{t) = it. Equation (|49)) turns into the parabolic 
cylinder equation. Its standard solution u{t, 0) has the 
asymptotic behavior: |w(oo, 0)p = exp(— 7rA^/2fi,e). Em- 
ploying it together with ([5^ and we arrive at the 
following number of molecules in the final state: 



N^{+oo) ^ F§N~^ [1 - exp (-7rAV2?ie)] 



(59) 



Eq. (|44|) implies that the maximal possible value of F^ is 
NsN/2. It corresponds to the complete transformation 
of atoms into molecules in the adiabatic regime e 0. 
Equation (|55)l looks exactly like the LZ transition prob- 
ability multiplied by an effective number of pairs. How- 
ever, in contrast to phenomenological thcoriesii^ and the 
perturbation theory result^, the coefficient in front of 
1 /e in the exponent does not depend on the initial atomic 
density. This theoretical prediction can be checked exper- 
imentally. The interaction and many-body effects influ- 
ence the effective number of pairs (pre-exponent), which 
is proportional to the square of the initial BCS conden- 
sate amplitude F{to). Finite temperature destroys a frac- 
tion of the initial Cooper pairs and decreases the molec- 
ular production. 

Finally, we consider an inverse process with no 
fermions, no BCS condensate and only the molecular 
condensate in the initial state: (&)(— oo) = \/N/2 and 
sweeping of the magnetic field in the opposite direc- 
tion. Then at the end of the sweeping the conden- 
sate density is determined by the LZ value: (6)(-|-oo) = 
y/ N/2 cxp(— 7rA^/2fi,e), whereas the absolute value of the 
BCS condensate amplitude {S-) can be found from the 
conservation law (|57p (for macroscopic condensate am- 
plitude we can neglect the non-commutativity of S+ and 

s-.y. 



\{S-)\' 



1 — exp 



ttA^ 
He 



(60) 



Notice the factor of 2 difference in the exponents of ([55]) 
and dSni)- 

An important problem is the reversibility of the pro- 
cess. Eq. (|49p for the amplitude b is time-reversible. 
It means that, if b{t) is a solution of ([49]) in a time in- 
terval (ti,t2), then b*{t) is also a solution of the same 
equation with e{t) replaced by e{t), where t denotes the 
time in the same interval passed in opposite direction. 
Unfortunately, it is not clear how to reverse the phase 
of the BEC in the experiment. If we simply start with 
the same number of molecules in the initial state, which 
we just before obtained after passing the FR, it does not 
guarantee that we reproduce in the end the state of the 
Fermi gas from which the molecules were obtained in 
the first half-cycle. It is obvious for the non-adiabatic 
process with not too large LZ parameter A^/e. In this 



9 



case to achieve reversibility we need not only to repro- 
duce correct ratio of number of molecules to the BCS 
condensate density, but also their relative phase which 
is large and quickly varying value. Even in the adia- 
batic situation the reversibility is scarcely reachable for 
two reasons. First, all levels except of the ground state 
have numerous crossings with other levels (see the previ- 
ous section) and the adiabatic approximation for them is 
invalid as it was indicated for more general situation by 
Polkovnikov and Gritse"\^. But even when we start with 
the ground state and sweep the magnetic field adiabati- 
cally, the asymptotic value of LZ process is reached only 
after very long time, since the amplitudes of "parasite" 
states decay slowly, as 1/t. The real experiment stops at 
some finite time, and corrections to the LZ result can be 
sufficiently large. They are unpredictable and represent 
a source of irreversibility in the experiment. 

VI. CONCLUSION 

First we estimate numerical values of different con- 
stants entering our theory from available experimental 
data. The numerical value g can be extracted from the 
experimental data on the magnetic field dependence of 
the scattering length a near the Fli^ using the well- 
known relation: g = h\J 47r(a — aQ)elm (ag is the scat- 
tering length far from resonance). On the other hand g 
can be estimated theoretically as (7 ~ e^,/ \/a^, where e/j / 
is the hyperfine energj*^. Both these estimates give for 
'*°K g ~ W>-'^^erg x cm^/^ and from cqs. (dH) and (fTT)) 
A ~ 3 X Wr'^K. However, eq. ([TT]) overestimates A by 
assuming that the limiting kinetic energy is A instead 
of it being much smaller. A more reliable estimate can 
be extracted from a comparison of eq. (j59p with experi- 
mental data by Regal et alJ^. The fitting gives the value 
A 10~^K for the broad resonance at Bq = 224.21G' in 
^''K. The cited measurements were performed at the finite 
temperature T ~ Tp/3,and therefore the corresponding 
value of A is underestimated in comparison to its zero- 
temperature value. Thus, a reasonable estimate for A is 
between lO^^K and lO^^K. In the cited experiment the 
magnetic field sweep amplitude was about 12G. It corre- 
sponds to an energy scale of about 10^'^ A', larger than 
A. 

In conclusion, we considered a cooled Fermi-gas in the 
magnetic field close the broad Feshbach resonance. We 
showed that for this strong coupling regime simplifying 
assumptions stemming from the broad resonance condi- 
tions allow one to solve both static and dynamic problems 
of the BEC-BCS conversion exactly. It was demonstrated 
that in this situation the single mode approximation is 
appropriate. The neglect of the fermion kinetic energy, 
being a rough approximation, gives qualitatively correct 
physical picture. With these two assumptions we derived 
the Global Spin Model Hamiltonian (|23p and its Hilbcrt 
space for the static case. For the dynamic conversion 
problem we solved the Landau-Zener problem for opera- 



tors. 

For the static problems we have found complete spec- 
trum, the number of molecules and its fluctuations, the 
amplitude of the BCS condensate and its correlation with 
the amplitude of the BEG condensate at different values 
of the detuning parameter e (see equations (|42|) . ()43p . 
(m, and (SSI)). 

One of the most surprising results is the big size of the 
Cooper pair in the momentum space and small its size in 
configurational space. Though we obtained this result in 
the framework of the GSM, it will remain in more pre- 
cise models since it is a consequence of a large extension 
of the hyperflne interaction in momentum space in com- 
parison with the Fermi momentum. Experimentally this 
theoretical prediction can be verified by the time-flight 
spectroscopy combined with the spin correlation analy- 
sis. After switching off the trap and the driving mag- 
netic field the Cooper pairs become unstable and decay 
forming the correlated pairs of particles with opposite ve- 
locities and parallel spins. The relaxation time is rather 
long since fermion pair collisions do not produce energy 
relaxation. The experimental estimate for the relaxation 
time is in the range from milliseconds to seconds. There- 
fore, it seems quite feasible to switch off the trap before 
the quasimolecules relax and observe the correlations of 
momenta and spins in runaway particles. 

For the dynamic problem we have considered molecule 
formation and dissociation when the magnetic field is 
swept across the broad Feshbach resonance. The re- 
sulting molecular production from initial fermions is de- 
scribed by LZ-like formulae with a strongly renormalized 
LZ gap independent of the initial fermion density. How- 
ever, the molecular production strongly depends on the 
initial value of magnetic field. Though being reversible 
in principle, the process is irreversible in practice since 
the preparation of the inverse process requires fixing of 
uncontrollable phases and slowly decaying corrections to 
the LZ asymptotics. Another experimentally verifiable 
prediction is the independence of the coefficient in front 
of 1/e in the LZ exponents for the molecular (|59|) and the 
BCS condensate productions of the initial density of 
atoms (molecules). 

Acknowledgments 

We acknowledge illuminating discussions with Victor 
Guraric, Leo Radzyhovsky, and A. Polkovnikov. Our 
thanks due to W. Saslow for attentive reading of the 
manuscript and constructive suggestions. Ar.A. grate- 
fully acknowledges financial support from Welch Founda- 
tion and NSF Grant PHY-0757992. V.P. is thankful to 
the Kavli Institute of Theoretical Physics, UCSB, for the 
hospitality extended to him during the Workshop on co- 
herence effects in cold gases and strongly correlated sys- 
tems. May - June 2007. This work was partly supported 
by the DOE under the grant DE-FG02-06ER46278 and 
KITP NSF Grant No. PHY05-51164. 



10 



^ J. Chwedenczuk et al., Phys. Rev. Lett. 93, 260403 (2004). 
^ K. Goral et al., Journal of Physics B: Atomic, Molecular 

and Optical Physics 37, 3457 (2004). 
^ F. H. Mies, E. Tiesinga, and P. S. Julienne, Phys. Rev. A 

61, 022721 (2000). 

* S. S. Inouye et al., Nature 392, 151 (1998). 

5 W. C. Stwalley, Phys. Rev. Lett. 37, 1628 (1976). 
^ E. Tiesinga, B. J. Verhaar, and H. T. C. Stoof, Phys. Rev. 
A 47, 4114 (1993). 

E. Timmermans, P. Tommasini, M. Hussein, and A. Ker- 
man. Physics Reports 315, 199 (1999), ISSN 0370-1573. 

* C. Chin et al., Phys. Rev. Lett. 90, 033201 (2003). 

^ J. CubizoUes et al., Phys. Rev. Lett. 91, 240401 (2003). 
1" E. A. Donley et al.. Nature 417, 529 (2002). 
" S. Durr et al., Phys. Rev. Lett. 92, 020406 (2004). 

J. Herbig et al.. Science 301, 1510 (2003). 

S. Jochim et al.. Science 302, 2101 (2003). 
" C. Regal et al.. Nature 424, 47 (2003). 

K. E. Strecker, G. B. Partridge, and R. G. Hulet, Phys. 

Rev. Lett. 91, 080406 (2003). 
1^ M. Greiner, C. A. Regal, and D. S. Jin, Nature 426, 537 

(2003). 

" E. Hodby et al., Phys. Rev. Lett. 94, 120402 (2005). 
M. W. Zwierlein et al., Phys. Rev. Lett. 91, 250401 (2003). 

1" N. R. Claussen et al., Phys. Rev. Lett. 89, 010401 (2002). 

2° D. M. Eagles, Phys. Rev. 186, 456 (1969). 

'^^ A. Leggett, in Modern trends in the theory of condensed 
matter (Springer- Verlag, Berlin, 1980), pp. 13-27. 
P. Nozieres and S. Schmitt-Rink, Journal of Low Temper- 
ature Physics 59, 195 (1985). 

Q. Chen, J. Stajic, S. Tan, and K. Levin, Physics Reports 
412, 1 (2005). 

C. A. R. Sa de Melo, M. Randeria, and J. R. Engelbrecht, 
Phys. Rev. Lett. 71, 3202 (1993). 

B. E. Dobrescu and V. L. Pokrovsky, Phys. Lett. A350, 
154 (2006). 



E. Timmermans et al., Phys. Lett. A 285, 228 (2001). 
^'^ J. Javanainen et al., Phys. Rev. Lett. 92, 200402 (2004). 
L Tikhonenkov, E. Pazy, and A. Vardi, Optics Communi- 
cations 264, 321 (2006). 

J. E. WiUiams et al, J. Phys. B 37, L351 (2004). 
^° V. Gurarie and L. Radzihovsky, Annals of Physics 322, 2 

(2007) . 

R. H. Dicke, Phys. Rev. 93, 99 (1954). 

D. Sun, A. Abanov, and V. L. Pokrovsky, Europhysics Let- 
ters 83, 16003 (2008). 

L. D. Landau and E. M. Lifshitz, Course of Theoret- 
ical Physics, Quantum Mechanics, vol. 3 (Butterworth- 
Heinenann, Linacre House, Jordan Hill, Oxford OX2 8DP, 
England, 2002), 3rd ed., translated from Russian by J. B. 
Sykes and J. S. Bell. 

E. M. Lifshitz and L. P. Pitaevskii, Course of Theoretical 
Physics, Statistical Physics, Part 2, vol. 9 (Butterworth- 
Heinenann, Linacre House, Jordan Hill, Oxford OX2 8DP, 
England, 2002), 1st ed., translated from Russian by J. B. 
Sykes and M. J. Kearsley. 

'^^ Y. Kagan and B. V. Svistunov, Sov. Phys. JETP 78, 187 
(1994). 

Y. Kagan, B. Svistunov, and G. Shlyapnikov, Sov. Phys. 
JETP 75, 387 (1992). 

C. W. Gardiner, P. ZoUer, R. J. BaUagh, and M. J. Davis, 
Phys. Rev. Lett. 79, 1793 (1997). 

C. W. Gardiner, M. D. Lee, R. J. Ballagh, M. J. Davis, 

and P. ZoUer, Phys. Rev. Lett. 81, 5266 (1998). 

P. W. Anderson, Phys. Rev. 112, 1900 (1958). 

A. Polkovnikov and V. Gritsev, Nature Physics 4, 477 

(2008) . 

C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 
92, 040403 (2004). 
■'^ For simplicity, we focus on a homogeneous system of vol- 
ume V 



